Mass inflation in a D dimensional Reissner-Nordstrom black hole: 
a hierarchy of particle accelerators ? 
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We study the geometry inside the event horizon of perturbed D dimensional Reissner-Nordstrom- 
(A)dS type black holes showing that, similarly to the four dimensional case, mass inflation also 
occurs for _D > 4. First, using the homogeneous approximation, we show that an increase of the 
number of spatial dimensions contributes to a steeper variation of the metric coefficients with the 
areal radius and that the phenomenon is insensitive to the cosmological constant in leading order. 
Then, using the code reported in [T] adapted to D dimensions, we perform fully non-linear numerical 
simulations. We perturb the black hole with a compact pulse adapting the pulse amplitude such that 
the relative variation of the black hole mass is the same in all dimensions, and determine how the 
black hole interior evolves under the perturbation. We qualitatively confirm that the phenomenon 
is similar to four dimensions as well as the behaviour observed in the homogeneous approximation. 
We speculate about the formation of black holes inside black holes triggered by mass inflation, and 
about possible consequences of this scenario. 



I. INTRODUCTION 



It is well established that the inner horizon of the 
Reissner-Nordstrom (RN) black hole is unstable and 
evolves into a curvature singularity when perturbed il}- 
I14j . The formation of this singularity is a consequence 
of a large relativistic counter-streaming produced by the 
perturbation in the vicinity of the inner horizon. Thus 
the black hole interior, when perturbed, effectively be- 
comes a particle accelerator, with centre of mass energy 
swiftly reaching Planckian regimes. One immediate con- 
sequence is that Planckian and trans-Planckian physics 
becomes relevant inside a black hole before the curvature 
scale of the unperturbed black hole becomes Planckian 
(in the latter case this should only occur close to the 
central curvature singularity of the classical solution). 
Trans-Planckian scattering is dominated by graviton ex- 
change [TS] and, for sufficiently small impact parameter, 
black holes should form [16|. Thus, the above picture of 
mass inflation leads us to the fascinating scenario where 
black holes could form inside black holes. 

In recent years, the physics of trans-Planckian scatter- 
ing has been intensively explored motivated by models 
of large P^THTO] or infinite [501 HI] extra dimensions. In 
these scenarios the fundamental Planck scale could be 
considerably smaller than the apparent four dimensional 
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Planck scale (10^^ GeV) and black holes could even po- 
tentially be formed in ongoing particle accelerator runs 
[2^ [25] . These black holes would be higher dimensional 
charged black holes, by virtue of the charge of the col- 
liding partons. Then, if these higher dimensional black 
holes also exhibit the same type of mass inflation instabil- 
ity as their four dimensional counterparts, we could have 
the amusing scenario of a (black hole) particle accelerator 
inside a (terrestrial) particle accelerator. 

In this paper we show that higher dimensional charged 
black holes exhibit mass inflation, quite analogous to 
their four dimensional counterparts (see |24j for a recent 
study of lower dimensional black holes) . More concretely 
we consider dynamical charged black hole solutions which 
generalise the standard Reissner-Nordstrom solution in 
D dimensional Einstein-Maxwell theory coupled to a cos- 
mological constant (Sections [Tl] - [rv] ). The role of the 
latter is shown to be quite irrelevant by considering the 
homogeneous approximation (Section |v|, which also al- 
lows the observation that the phenomenon is faster (in 
terms of the areal radius) as one increases the dimension, 
for the same (perturbation) energy density placed at the 
same areal radius. We then perform full numerical sim- 
ulations of a compact wave packet falling into the black 
hole (Section VI), using the code described in [T]. To 
make a comparison between the different dimensions, we 
perform runs in D = 4, 5, 6, 7 varying the amplitude of 
the pulse in such a way that the relative black hole mass 
variation is constant, approximately 25%. We exhibit 
the dynamics of the inner (and also of the outer) hori- 
zon, which is qualitatively quite similar in all dimensions. 
We also display the behaviour of the mass function and 
areal radius inside the horizon, from which a qualitative 
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agreement with the observations made for the homoge- 
neous approximation may be inferred. 



II. ACTION AND FIELD EQUATIONS 

Consider the Einstein-Maxwell action with a cosmo- 
logical constant and minimally coupled to a scalar field 
in a D dimensional spacetime: 



S = J d^xy/-g (^R 



{D-l){D-2) 



A + 167r£ 



M 



(1) 

where R is the Ricci scalar and A is the cosmological 
constant. Lm is the matter Lagrangian, which we will 
assume to be the sum of the usual Maxwell contribu- 
tion, Lp, and the contribution from a self-gravitating 
real massless scalar field ip: 
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where = FaisF"^ and Fap is the Maxwell tensor. Our 
signature choice is — h -I- -t- . . . . The energy-momentum 
tensor of the matter fields is given by ^T^,y + '^T^,^, where 



Ft - — ( F F " - -a F^ \ 
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(3) 
(4) 



The equations of motion derived from ([T]) with ^ are 
the gravitational equations 

G,, + ' ^^^^ ' ^h g^, = 8^(^r^, + ^T^,) ,(5) 

the matter equations 

□(^ = , d*F = , (6) 
where ★ denotes Hodge dual. 



III. D DIMENSIONAL 
REISSNER-NORDSTROM-(A)DS 

A well known solution to this system is that of a _D di- 
mensional spherically symmetric (i.e. with spatial isom- 
etry SO{D — 1)) charged black hole in asymptotically 
(Anti)-de Sitter spacetime, first discussed by Tangherlini 
[ig (see also [2B]). The fields read 
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where dQo- 
and 



is the line element on the {D — 2)-sphere 
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We have chosen the coefficient of the cosmological con- 
stant term in the gravitational action such that the A 
term in the last equation is dimension independent. The 
two metric parameters M and q are related to the ADM 
mass and charge by 
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wheve Ad-2 is the area of a (Z? — 2)-sphere of unit radius. 

The solution has Killing horizons of the d/dt Killing 
vector field at radial distances r obeying 
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If A = 0, the solutions are a simple generalisation of the 
four dimensional case; the horizons are located at 



M ± ^AP - g2 



(13) 



Thus, the higher dimensional Reissner-Nordstrom- 
Tangherlini black hole still has an inner, at r = r_, and 
an outer horizon, at r = r+, as long as M > \q\. In this 
paper we shall always take M = 1 and q — 0.95 for the 
unperturbed black hole. 

If A ^ one cannot, in general, provide an analytical 
solution for (12 1. Qualitatively, if A < 0, the spacetime 



is asymptotically Anti-de-Sitter and there are no further 
Killing horizons; if A > 0, the spacetime is asymptotically 
de-Sitter and there may be further Killing horizons of 
cosmological nature. 



IV. ANSATZ WITH SPHERICAL SYMMETRY 
IN DOUBLE NULL COORDINATES 

To solve numerically the coupled system ^ and ([6| 
is, in general, quite a difficult task. In fact, numeri- 
cal relativity in higher dimensions is a very recent field 
(see eg. [27H33] for general formulations with applica- 
tions). Here we are interested in spherically symmetric 
black holes, perturbed in such a way that this symmetry 
is preserved. This makes the system tractable. Thus, we 
take a spherically symmetric ansatz written in double- 
null coordinates: 



ds^ = -2e2'"("'")rfudt) + P{u, v)dnD- 
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where u and v are taken to be ingoing and outgoing re- 
spectively. For the unperturbed sohition of the previous 
section, <y9 = 0, e^'"^"'") = /(r) with 

u=^{t-r*), „=.-L(< + r*), (16) 

and dr/dr* — f{r). 

The Maxweh equations ^ are simply solved to yield 

FuvT^^^e^^" = constant = Q , (17) 



where 



2A 



D-2 



(18) 



is the electric charge. The electric field is therefore purely 
radial, as expected from spherical symmetry. The scalar 
field equation (|6| gives 
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The uu, vv, uv and transverse components of the Ein- 
stein equations ([s]), which are the only distinct non-trivial 
ones, give rise to the following equations (the left hand 
sides of these equations are components of the Einstein 
tensor G^i,): 
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V. HOMOGENEOUS APPROXIMATION 

In the homogeneous approximation one considers only 
r dependence in the line element (see jl] for a discussion 
of this approximation and references). That is we keep 
the symmetry with respect to the d/dt Killing vector, 
which is spacelike in between the inner and outer horizons 
of the unperturbed solution. The spherically symmetric 
homogeneous line element is given by 



ds"^ = gu{r)dt^ + grr{r)dr'^ + r'^dVt 
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The black hole charge will produce a purely radial elec- 
trostatic field. The Maxwell equations then yield 
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The scalar field admits solely a radial dependence, and 
the Klein-Gordon equation may be written as a first order 
equation 
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where the prime denotes radial derivative. The tt and rr 
components of the Einstein equations read, respectively 
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A. Analytic analysis for the homogeneous 
~ approximation 

(23) The scalar field we have been considering may be re- 
garded as a perfect fiuid (see eg. [SJ) with 



For D = 4, A = 0, these equations reduce to the ones in 
[T] with constant Brans-Dicke scalar. 

Equations (22)-(23) are the evolution equations, that 



determine the dynamical variables in the future of the 
hypersurfaces where initial conditions are set (see Fig. 1 
in [1]). In equations (20)-(21), the contribution from the 



energy-momentum of the scalar field is 
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respectively, which represent, physically, the flux of the 
scalar field through surfaces of constant v and u, i.e. out- 
flux and influx. We shall impose that, initially, only in- 
flux exists. However, outflux is inevitably produced by 
scattering off the spacetime geometry. 



p = p = AnC 
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(30) 



In the homogeneous approximation = '/'('')■ Then, 
energy momentum conservation for the scalar field yields 
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where pi, gtu and are all integration constants; pi is 
the density at the surface r — ri, where gtt is gtti- In the 
mass infiation region, we expect the variations of gtt and 
grr to be fast. Thus, we approximate (28) and (29) by, 
respectively 
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Observe that the cosmological constant becomes irrele- 
vant, in leading order. Moreover, we conclude that, as in 
four dimensions, the radial and temporal components of 
the metric should be approximately proportional 
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It follows from this result, together with (31) and (33) 
that 
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component is more abrupt as the spacetime dimension 
increases. Fixing and pi, d\ngrr/dr oc {D — 3)^{D — 2) 
if ri — 1. This can be confirmed in Fig. [2j where we plot 
the radial metric component for two different values of pi 
and for D — 5,6, 7. Indeed Fig. [2] shows that the slopes 
increase with D, for a fixed pi. Note that the radial co- 
ordinate is the areal radius, therefore having an invariant 
geometric meaning. Thus it is the appropriate variable 
to compare the mass inflation phenomenon in different 
dimensions. In order to make this comparison, however, 
we fixed = 1 in all dimensions. In the next section, we 
shall fix the same physical observable, namely the same 
relative mass variation for all dimensions. 



where 
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is roughly constant in the mass inflation region. The 
transition from the RN phase (where g^- and gtt given 
approximately by Eqs. ([7| and ([9| with —grrgtt = 1) 



to the mass inflation phase (where Eq. (34) is approxi- 
mately valid) is rather sharp and can be defined by the 
equality of the two terms on the r.h.s. of Eqs. (28) 
and ( 29 ) . Taking into account that at the transition 
{—grr/ guY^^ — —grr — ^/ gtt, one finds that, in the mass 
inflation region. 
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Taking the logarithm of g^r and expanding it up to linear 
order in r — r_ in the mass inflation region (r < r_ ) one 
gets 
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FIG. 1: Variation of the radial and temporal metric compo- 
nents with the areal radius normalised by the value at the 
inner horizon, for a D = 5 charged black hole. 
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Using Eqs. (jllj, ([13f, (jSG)), (|37|) and ([SSf one finally 
obtains 
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This leads to the near straight lines in the mass infla- 
tion region in Fig. [T] where we plot the variation of the 
radial and temporal metric components for D = 5. As 
in the four dimensional case. Fig. [T] shows that mass 
inflation is more abrupt for smaller values of pi. This 
happens because dhigrr/dr oc \/pi in the mass infla- 
tion region which results in a larger slope of In grr {t) for 
smaller values of pi. In higher dimensions the behaviour 
is similar but the relative variation of the radial metric 
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FIG. 2: Comparison of the radial metric component behaviour 
for D — 5, 6, 7 and for two different values of pi, in terms of 
the areal radius normalised by the value at the inner horizon. 
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VI. FULL NUMERICAL SIMULATIONS 

Horizons 



0.00 
5 



D=i, A =0.275 
D=5, A =0.125 
D=6, A =0.065 
D=7, A =0.0343 
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FIG. 3: Evolution of the black hole mass as it interacts with 
the pulse. The amplitude was chosen, in each dimension, so 
that the final black hole mass variation is approximately 25% 
for all cases. 



We shall now present the results we have obtained us- 
ing the numerical code reported in [1] generahsed to D 
dimensions and the compact pulse used therein, which is, 
along the initial null segment u = uq^ 



f,viuoi v) = A sin TT 



Av 



vo < V < vq + Av , 



(41) 

and vanishing outside the interval vq to vq + Av, with 
an interval Av — 1. With this code we have performed 
fully non-linear numerical simulations in I? = 4, 5, 6, 7. 
In order to make a meaningful comparison between the 
different dimensions, our strategy was to fix the same rel- 
ative mass variation, AM/M, after the interaction with 
the pulse, for the different dimensions. In order to achieve 
this we had to vary the amplitude of the pulse A by an 
iterative procedure. In Fig. [3] we show the evolution 
of the mass of the black hole for the values of the am- 
plitude that were found to give, for each dimension, a 
relative mass variation of about 25%. 

In order to follow the black hole mass, we computed 
the (outermost) apparent horizon at each null slice, i.e. 
the set of points {ua,Va) such that 



dr 
dv 



iUa,Va) = 



(42) 



In Fig. |4]we exhibit the apparent horizon evolution (both 
inner and outer). Then, the black hole mass was given 
by the Misner-Sharpe |35j mass function evaluated at this 
apparent horizon. The Misner-Sharpe mass is the total 
effective mass inside a sphere of radius r{u,v), and is 
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FIG. 4: Evolution of the black hole apparent horizons. For 
each dimension the outer (top curve) and inner (bottom 
curve) are plotted. 
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In Fig. [5] we have plotted the Misner-Sharpe mass 
function in the mass inflation region. The plots show 
that larger masses occur for lower values of advanced and 
retarded time, as the dimension increases. Although this 
is a gauge dependent statement, we emphasise that both 
t and r that define u and v at infinity are comparable 
coordinates in different dimensions (proper time of the 
observer at infinity and areal radius). 

It is tentative to conclude that the smaller time scale 
for the mass variations (and also curvature as we have 
checked) observed in our fully non-linear numerical sim- 
ulations, for increasing dimension, is consistent with the 
more abrupt radial variation of the metric coefficients, 
observed in the homogeneous approximation as the di- 
mension increases. This is qualitatively confirmed by a 
simultaneous analysis of Fig. [5] and Fig. [6] where the 
areal radius is plotted in the same (u, v) space. Con- 
sider, for instance, a v — constant line (say v = 9), for 
D = 7. From Fig. [5] it is seen that a variation of 10 
orders of magnitude of the mass occurs in an interval of 
Au ~ 0.1 — 0.2 for u ~ 15.7. Comparing with Fig. [6] 
we see that the variation of the areal radius in the same 
interval is around Ar ~ 0.1 — 0.2. In this r interval, there- 
fore there is a mass variation of 10 orders of magnitude. 
A similar analysis, along the same v = constant line in 
a lower D shows a considerable smaller mass variation 
for the same r interval. For instance, in D = 4 the mass 
variation is of 1 — 2 orders of magnitude in the same r 
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FIG. 5: Misner-Sharpe mass function for a perturbed RN black hole in D = 4, 5, 6, 7, with the amplitude of the perturbation 
chosen such that the relative mass variation of the black hole is the same in all dimensions. As the dimension increases larger 
masses are achieved for lower w, v. 



interval. 

One must be careful in making this connection, how- 
ever, since both (m, v) and r are simply coordinates and 
therefore gauge dependent; moreover in each case one 
is fixing different quantities. In particular, in the ho- 
mogeneous approximation one is fixing the initial energy 
density at a given point, rather than the total energy in 
the perturbation. 



VII. CONCLUSIONS 

If micro black holes form by trans-Planckian scattering 
collisions in the present or in a future generation of par- 
ticle accelerators, they should undergo, after formation, 
a balding phase, in which they loose gravitational (and 
other interactions') multipoles and approach a stationary 
solution, a semi-classical phase, in which the black hole 
decays by Hawking radiation, with concrete and well de- 
fined signatures that could be observed in particle exper- 



iments |36j , and finally a mysterious (with present day 
physics) Planckian phase [37] . 

For the scenarios proposed in [T7H2T] , these black holes 
are brane world black holes. But if their mass is suffi- 
ciently above the fundamental Planck scale and their size 
is sufficiently smaller than the scale of the large extra di- 
mensions, then, in the absence of charge, they should be 
well approximated by the classical solution for uncharged 
black holes in infinite D dimensions. The charged case 
is more involved, since the electromagnetic field should 
propagate only on the brane (which must have a length 
scale smaller than the TeV) and therefore in 3 -I- 1 rather 
than in D dimensions. Modelling such black holes as a 
D dimensional RN black hole is the simplest toy model, 
which we have explored in this paper. 

According to the result of this paper, the brane world 
black holes will undergo mass inflation and, potentially, 
form new smaller black holes inside of themselves, at 
least if the classical/semi-classical stage has a sufficiently 
long time scale. One could think this would lead to a 
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D=4, ^=0.275 



D=5, yl=0.125 





23=6, ^=0.065 



D=7, ^=0.0343 





FIG. 6: The radial coordinate r as a function of the nuU coordinates u and v close to the inner ingoing apparent horizon for 
-D = 4, 5, 6, 7 perturbed RN black hole with M — 1 and q = 0.95 and the amplitudes that guarantee a total mass variation of 
AM/M ~ 0.25 . Contours are plotted only for radii in the range of the color bar. 



sort of fractal, self-similar structure of black hole gener- 
ations. Remember, however, that the minimum mass for 
the black holes is the (fundamental) Planck mass, which 
therefore should limit the number of generations to one. 

One other possible scenario is as follows. These ex- 
tremely energetic regimes are also expected to be at- 
tained in the early universe. Although the relevant 
physics is still not fully developed, it has been postu- 
lated that an inflationary universe might naturally arise 
in such extreme conditions, solving some of the problems 
of the standard cosmological model. Hence, the possibil- 
ity that mass inflation inside black holes may also trigger 
the creation of new expanding universes should be consid- 
ered. In [55] it was assumed that quantum effects could 
remove the black hole singularities and be responsible for 
the birth of new cosmologies with slightly different values 
of the cosmological parameters, thus leading to a selec- 



tion criteria based on the maximization of the number of 
black holes. The fact that Planckian and trans-Planckian 
regimes are attained in an extended mass inflation region 
much before the central singularity and the role played 
by the number D of spacetime dimensions on the corre- 
sponding dynamics must therefore be taken into account 
in such studies. 
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